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Both CDF and D0 experiments have been searching for evidence of physics beyond the stan- 
dard model (SM) using the Tevatron pp collider at Fermilab. We report on recent searches in 
the B flavor sector, especially decays via flavor changing neutral current processes (FCNC), 
_B°^j e+p" and -B^ M"*^/^") at the Tevatron. 



1 Introduction 

The branching fraction of a rare decay mode is an interesting quantity to measure because the 
contribution from physics beyond the SM may be sizable in rare decay modes. In order to be 
able to observe rare heavy flavor decays it is essential to produce a sufficient number of bottom 
hadrons. The bottom anti-bottom production cross-section a{bb) at the Tevatron is O(IO^) 
larger than production in e'^e" colliders at the T(4s) or energy scale. The large production 
of all kinds of 6-hadrons at the Tevatron offers the opportunity to study rare decays also in 
the -Bg and 6-baryon sectors. On the other hand the inelastic cross section is 10^ times higher 
than (t(66) requiring very selective and efficient triggers. Therefore, interesting events must be 
extracted from a high track multiplicity environment, and detectors need to have very good 
tracking, vertex resolution, wide acceptance, good particle identification and highly selective 
triggers. 

In this article we report on searches for rare B decays via FCNC process, B^^^^ — > and 
B^ — > ^'^ jjL~ , using the Tevatron pp collider at Fermilab. 



2 ^l*^^ e+/x- Decays 

In the SM, lepton number and lepton flavor are conserved; therefore the decays such as B^^^^ — > 
e^/i~ are forbidden. However the observation of neutrino oscillations indicates that lepton flavor 
is not conserved while lepton flavor violating (LFV) decays in the charged sector have not been 
observed yet. The grand-unification theory by J. Pati and A. Salam predicts a new interaction to 
mediate transitions between leptons and quarks via exchange of spin-1 gauge bosons, which are 
called Pati-Salam leptoquarks (LQ), that carry both color and lepton quantum numbers^. The 
lepton and quark components of the leptoquarks are not necessarily from the same generation 
and the decays B^ e'^fi" and B^ — > e~^fi~ can be mediated by different types of leptoquarks. 

CDF reported on searches for the LFV decays B^_,^ — > e~^fj,~ using 2 fb~^ of Run II data^. 
Data sample used in the searches was taken with the two displaced-track trigger, in which two 



oppositely-charged tracks are required to have a transverse momentum px > 2 GeV/c and an 
impact parameter^ 0.1 < do < 1 mm, and it is also required that the scalar sum of the transverse 
momenta of the two tracks is greater than 5.5 GeV/c, the difference in the azimuthal angles of 
the tracks is 20° < Ay? < 135°, and a transverse decay length^ is Lxy > 200 //m. In the off-line 
analysis, additionally the B^^-^ isolation, the pointing angle and the transverse decay length were 
required to be consistent with those of B^^^ — > e"*"/!" decays. These thresholds were optimized 
in an unbiased way to obtain the best sensitivity for the searches ^. To identify electrons, both 
the specific ionization (dE / dx) ^ measured in the central drift chamber and the transverse and 
longitudinal shower shapes as measured in the central electromagnetic calorimeter were used. 
The electron identification efficiency is ~ 70% while the muon identification is fully efficient with 
> 2 GeV/c in the central muon detector. Search windows in e^/i~ invariant mass were defined 
to be (5.262-5.477) GeV/c^ for B^ e+^i' and (5.171-5.387) GeV/c^ for B° e+^i'. These 
correspond to a window around the nominal values of the B^ and B^ masses^ of approximately 
=b3(Tm. The background contributions considered include combinations of random track pairs 
and partial B decays that accidentally meet the selection requirement and hadronic two-body B 
decays in which both final particles are misidentified as leptons. The random track contribution 
was evaluated by extrapolating the normalized number of events found in the sidebands to the 
signal region. The double-lepton misidentification rate was determined by applying electron 
and muon misidentification probabilities to the number of two-body decays found in the search 
window. Figure 1 shows the invariant mass distribution for e^fi~ candidates. One event in the 
B'^ mass window, and two events in the B^ mass window were observed. These numbers are 
consistent with the estimated total background of 0.8 it 0.6 events in the B^ search window, and 
0.9 ± 0.6 in the B^ window. B^ — > K^Tr~ decays were used as a reference channel to set a limit 
on B{B°^^s^ e+fi-). Using the same selection criteria except lepton identification, 6387 it 214 
events of B^ K~^-k~ decays were observed as shown in Fig. 1. Resulting upper limits were 
B{B^ e+fi-) < 2.0 (2.6) xlO"'^ and 23(5° ^ e+^u") < 6.4 (7.9) xlO"*^ at 90% (95%) confidence 
level (C.L.). Using the limits on the branching fractions, the masses of the corresponding 
Pati-Salam leptoquarks were calculated to be Mlq(S° ^ e+n') > 47.8 (44.9) TeV/c^ and 
Mlq(S° e^fi') > 59.3 (56.3) TeV/c^ at 90 (95)% C.L. These are the best limits in the world 
to date. 
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Figure 1: ^ invariant mass distribution(left) and the reference channel B^' ~* K^-r (right) in data. 



3 S° ^ Decays 

Branching fraction of the pure leptonic FCNC processes like B^ —>■ fJ,~^fJ,~ in the SM is heavily 
suppressed B{Bg IJ'~^fJ'~) = (3.35±0.32) x 10~^. B^ — > decay is further suppressed by 

a factor |Vj^/VtsP in the CKM matrix elements leading to a SM predicted branching fraction of 
0(10-1°). The decay amplitude of 5° — > fi^fJ- can be enhanced significantly in some extensions 



of the SM. For instance, in the type- II two-Higgs-doublct-modcl (2HDM)-'-^ the branching fraction 
is proportional to tan^/3, where tan/3 is the ratio between the vacuum expectation values of the 
two neutral Higgs fields. In the minimal super-symmetric SM (MSSM) the dependence on tan/3 
is even stronger, oc tan^/3. Current observation of this decay at the Tevatron would necessarily 
imply new physics. 

Current world best limit on B{B^ is given by CDF using 2 fb^^ of Run II data^'^. 

The sensitivity of the analysis was improved significantly from the previous result by increasing 
the integrated luminosity of the event sample, using an enhanced muon selection, employing a 
neural network (NN) classifier to separate signal from background, and performing the search in 
a two dimensional grid in dimuon mass and NN space. The observed event rates were consistent 
with SM background expectations. Extracted 90% (95%) C.L. limit was B{Bg — > < 
4.7 (5.8) X 10-^ 

D0 reported a new expected upper limit on B{Bg M^M~) using approximately 5 fb~^ of 
dimuon trigger data^^. The data sample was split into three subsamples. Run Ha, Run Ilb-I and 
Run nb-II, based on the data taking period. Roughly the integrated luminosity of Run Ha data 
is 1.3 fb~^. Run Ilb-I is 1.9 fb~^ and Run Ilb-II is 1.6 fb~^. The three subsamples were treated as 
three different and independent analyses, but the final upper limit was combined from the sepa- 
rate analyses. To separate signal from background a boosted decision tree (BDT) classifier was 
constructed, where five input variables of the Eg meson (isolation, transverse momentum, trans- 
verse decay length significance, impact parameter significance and logarithm of vertex prob- 
ability) were used. Figure 2 shows the dimuon invariant mass distributions after applying the 
BDT cut, where only events not in the blinded region ±3o-„ (am = 0.115 GeV/c^) around the B° 
mass are shown. The same BDT cut was applied on the normalization channel, — > J/ipK'^ , 
and the following numbers of events were found: 1847±49(stat.) it 115{syst.) events in Run Ila 
data, 2188 ± 52{stat.) ± 123{syst.) events in Run Hb-I data and 1683 ± 46{stat.) ± 112{syst.) 
events in Run Ilb-II data. The systematic uncertainty is coming from our parametrization of 
the mass shape. The random track contribution in the signal region ±2.50"^ was estimated 
by extrapolating events in the sideband region to the signal region. In addition, possible non- 
negligible contributions of misidentified B^ K^K~ and B^ — > K'^-k" were estimated. The 
total background events in the search window in Run Ha, Riin Ilb-I, and Run Ilb-II data sets 
were estimated to be 2.16 it 0.62, 3.73 ± 1.07 and 2.15 ± 0.63, respectively. The expected SM 
yields of 5° l^^l^-^ events in Run Ha, Run Ilb-I, and Run Ilb-II data sets are 0.192 ± 0.034, 
0.193 lb 0.034 and 0.139 it 0.025, respectively. Aside from the background uncertainty, the largest 
uncertainty common to the three data sets, 15.2%, comes from the fragmentation ratio between 
B^ and Bg. Assuming no signal counts (background only) in the signal region, an expected 
upper limit on the branching fraction at the 90%(95%) C.L was computed. The number of 
observed events was set to the number of background events, 2 events for Run Ha, 4 events for 
Run nb-I and 2 events for Run Ilb-II. In this calculation, it was assumed that there are no 
contributions from B^ — >■ ji^ decays, where the decay is suppressed by jV^w/Vt^p « 0.04. The 
expected 90% (95%) upper limits for the branching fraction were 7.6 (9.4) x 10~^ for Run Ila 
data, 9.9 (11) x 10"^ for Run Hb-I data and 10 (13) x 10"^ for Run nb-H data. The combined 
upper limit is then 4.3 (5.3) x 10~^ at the 90%(95%) C.L. This sensitivity is comparable with 
the best published upper limit from CDF^^, and improves the previous D0 result by a factor 
of two. Further work to understand and reduce the background as well as to include more data 
is ongoing before opening the search box. 

4 Conclusion 

CDF has reported new upper limits on B{B^g-^ e+/x^) using 2 fb"-*^ of data, and D0 has 
reported an new expected upper limit on B{B^ —>■ iJL^iJt~) using approximately 5 fb~^ of data. 
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Figure 2: Dimuon invariant mass distributions. Search box remains blinded. 

With the continuously increasing amount of data provided by the Tevatron and improvements 
of the analyses measurements of rare B decays provide new insight into the properties of the 
FCNC decays, which allows for improved tests of the SM that could guide us to new physics 
scenarios. 
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